Energy-dependent effective interactions for dilute many-body systems 
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We address the issue of determining an effective two-body interaction for mean-field calculations 
of energies of many-body systems. We show that the effective interaction is proportional to the 
phase shift, and demonstrate this result in the quasiclassical approximation when there is a trap- 
ping potential in addition to the short-range interaction between a pair of particles. We calculate 
numerically energy levels for the case of an interaction with a short-range square- well and a harmonic 
trapping potential and show that the numerical results agree well with the analytical expression. 
We derive a generalized Gross-Pitaevskii equation which includes effective range corrections and 
discuss the form of the electron-atom effective interaction to be used in calculations of Rydberg 
atoms and molecules. 

PACS numbers: 03.75.-b, 34.50.-s, 34.60. +z 
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I. INTRODUCTION 

Ultracold gases represent an ideal environment to 
study fundamental processes such as interparticle colli- 
sions and molecule formation. These systems are gener- 
ally dilute, in the sense that the mean distance between 



particles 



-1/3 



is typically larger than the range R of 



the interactions: under these conditions many-body en- 
counters are rare, and interactions can be satisfactorily 
modeled by two-body collisions. 

Atomic interaction potentials have a complicated 
structure and are generally not known exactly but, as 
shown by Fermi [lj, the long- wavelength, low-energy 
properties of the two-body system can be reproduced 
exactly if one replaces the true potential by a suitable 
boundary condition for the relative wave function at the 
origin, 



ip(r) oc 1 



(1) 



where r is the relative coordinate vector. This boundary 
condition depends on a single parameter, the scattering 
length a = — lim^o S/k, where 5 is the s-wave scattering 
phase shift. When employed in the two-body problem, 
this leads to an energy shift given by 



2nh 2 a, , . , 

AE= |V(0)| 2 , 

M 



(2) 



where ^>(r) is the relative wave function in the absence of 
two-body interactions and [i is the reduced mass of the 
two particles. This result leads to the interaction term 
in the Gross-Pitaevskii (GP) equation for the mean-field 
wave function of a Bose-Einstein condensed gas. The 
standard Fermi treatment is justified as long as the rel- 
ative momentum k is so low that k\a\ <C 1, but to deal 



with phenomena at higher energies the theory must be 
improved. 

An extension of the Fermi pseudopotential was pro- 
posed in Refs. 0, 01- where the authors introduced a 
more general contact pseudopotential which depends on 
energy and which includes all partial waves I. Its s-wave 
component is given by 



Tr , 2nh 2 tan 5 , drijj 



(3) 
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The solution of the Schrodinger equation for the pseu- 
dopotential agrees with that for the actual potential in 
the region where the actual potential vanishes. The 
magnitude of the pseudopotential is specified in terms 
of the tangent of the phase shift, tan<5(fc), where the 
wave number k must be taken to be that in the ab- 
sence of the potential. As shown in 0, 0], when the 
Schrodinger equation is solved for this pseudopotential 
one obtains eigenenergies and wave functions that repro- 
duce very accurately the results obtained from integrat- 
ing the Schrodinger equation directly for the actual mi- 
croscopic potential. Following the suggestion of Ref. @], 
such a pseudopotential has been employed in studies of 
Rydberg atoms and molecules 0, EJJ- However, when 
used in mean-field calculations of the energy, it yields di- 
vergent energy shifts when a phase shift becomes close to 
7r/2 (modulo 7r). A further proposal for the relationship 
between an effective interaction and the phase shift has 
been made in the context of deriving a generalization of 
the GP equation, where it has been argued that the en- 
ergy shift should be proportional to the real part of the 
forward scattering amplitude, i.e. AE oc sin S cos 5 [ill ]. 

In this paper we explore the relationship between en- 
ergy shifts and phase shifts. We shall argue that the 
generalization of Eq. <|2j to nonzero k is to replace the 
scattering length by —5/k, a result demonstrated long 
ago for particles in the absence of a trapping potential 
|12 |. In Sec. II we first demonstrate this result for two 
particles whose relative motion is confined to lie within a 
sphere, and then show that in the quasiclassical approx- 
imation it also holds if there is an additional trapping 



2 



potential that depends on the relative coordinate. Fol- 
lowing that we calculate numerically the energy shift for 
an interaction consisting of a short-range square well plus 
a harmonic trapping potential and show that the analyt- 
ical formula fits the data well. In Sec. IIII1 we derive a 
generalization of the GP equation that includes the ef- 
fective range of the interaction through a derivative term 
and compare results with predictions based on the pre- 
scription of Ref. [JjJ- In Sec. IIVI we analyze effective 
interactions for ultracold Rydberg atoms and molecules. 
We summarize our results in Sec. 



II. ENERGY SHIFT 

The relationship between energy shifts and phase shifts 
for a particle interacting with a static impurity was an- 
alyzed long ago (see e.g., |l2j), and we briefly review 
the argument in the context of the two-particle problem. 
Consider two particles, with reduced mass fx, interacting 
via a spherically-symmetric potential. In the absence of 
interaction, the relative wave function for an s-state is of 
the form 



A 



sin(fcor) 



(4) 



For definiteness, we imagine the relative motion to be 
confined by a sphere of radius R, and we impose the 
boundary condition that the wave function must vanish 
at r = R. This implies that the wave number in the 
absence of interaction is k$R = nir, and normalization 
of the wave function gives A = (21TR)- 1 / 2 . In the pres- 
ence of a short-ranged interaction V sr (r), which we shall 
assume vanishes more rapidly than r _1 for large r, the 
asymptotic wave function will have the same form with 
a phase shift: 



V>(r) = A 



sm(kr + S) 



(5) 



To satisfy the boundary condition at r = R, the wave 
number must now obey the equation kR + S = nir, which 
implies that the wave vector is shifted by an amount 
Afc = k — ko = —6/R. The energy shift is then given 

by 



AE : 



-k Q Ak 



(6) 



Thus the energy shift due to interparticle interaction is 
proportional to the phase shift S. In the limit of zero 
energy, scattering theory (see e.g., shows that the 
s-wave phase shift is proportional to the wave vector 
S ~ —ka and one recovers the well-known result that the 
effective interaction has a contact form, with strength 
Uq = 2nti 2 a/[i. 



A. Presence of an external confining potential 

One may ask whether the result © applies in the pres- 
ence of a trapping potential. We therefore consider the 
same problem as above, but with an additional external 
potential V cx (r) for the relative motion We shall 

assume that V cx (r) increases with increasing r, and we 
shall impose the boundary condition that the wave func- 
tion tends to zero for large r. 

The basic assumption we shall make is that the trap- 
ping potential V cx (r) varies negligibly over both the range 
L of the two-body potential V sr (r) and over length scales 
~ \S/k\, which we shall show will play the role of an 
energy-dependent scattering length. In addition, we as- 
sume that Vcx(r) varies sufficiently slowly in space and 
that the energy of the state is sufficiently high that we 
may employ the quasiclassical approximation. As usual, 
it is convenient to work with the quantity \ = i^Pi an d 
in terms of it, the Schrodinger equation becomes 



2 M dr 2 X 



(7) 



and we may take x to be real. The quasiclassical approx- 
imation for the potential V cx (r) is 



A 



•yw) 



dr' \ 



1 2 t i[E - V cx (r')} 



. (8) 



where p(r) — y/2p,[E — V ex (r)] is the classical relative 
momentum of the two particles. The normalization con- 
stant A is fixed by requiring that the norm of ip, i.e. the 
volume integral of \ip\ 2 out to the classical turning point 
r = r c where V cx (r c ) = E, be unity. Since in Eq. ||HJ 
p^ 1 / 2 changes slowly over a period of oscillation of the 
sine function, we obtain 



dr i> 2 {r) = 2-kA 2 



dr 
p(r) 



1. 



(9) 



By choosing the zero of energy such that T4 X (0) = 0, the 
amplitude of the wave function when the positions of the 
two particles coincide is found to be 



^(0) 



2/iE/h 4 



^2tt J rc drp 1 (r) 



(10) 



The phase of the semi-classical wave function obeys 
Bohr's quantization rule 0] 



dr 



l2fi[E-V m (r)] 



= (n + a)n, 



(11) 



where n is a positive integer and a is a constant that 
depends on the nature of the potential in the vicinity of 
the classical turning point. The presence of the potential 



3 



V sr (r) induces an energy shift AE, and asymptotically 
the wave function acquires a phase shift S that satisfies 



dr 



l2n[E + AE~V cx {r)} 
h 2 



-8(E + AE) = (n + a)ir, 
(12) 

where r' c is the appropriate classical turning point in the 
presence of the short-range interaction. By taking the 
difference between Eqs. 11211 and llllt . expanding the in- 
tegral to first order in AE, and making use of Eq. (1101 . 
one obtains the result 



AE=- 2 -^ 6 M±^\m\\ 

2jjLE/fi 2 



(13) 



in agreement with Eq. 10. The difference between r c 
and r' c plays no role since the integrand vanishes at the 
turning point. 



B. A numerical example 

We now perform numerical calculations of the energy 
shift for a model potential. For the trapping poten- 
tial, we choose a harmonic potential of frequency ui, 
V cx (r) — fiuj 2 r 2 /2, while for the particle-particle inter- 
action we take an attractive spherically-symmetric po- 
tential which is equal to a constant, Vq < 0, for r < L 
and zero otherwise. The problem is similar to that con- 
sidered in Refs. 0,IHI3> but with a simplified short-range 
interaction. For r < L, the solution with energy E that is 
regular at the origin is given by Eq. Q with wave vector 

k = ^2^(E-V )/h 2 , 



i>in(r) oc 



sin(fcr) 



(14) 



For r > L, the wave function is the general solution of 
the Schrodinger equation for the harmonic potential, 
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Figure 1: Energy levels for the model potential described in 
the text as a function of phase shift, which is given by Eq. 
Ultifl . Levels are indexed by the principal quantum number n. 
The symbols indicate the results of numerically solving the 
matching condition for the Schrodinger equation at r = L, 
The values of the energies E(V BI ) = E(0)+AE(V BI ) calculated 
using Eq. 11311 (continuous lines) are indistinguishable from 
the values (symbols). As a comparison, we also show the 
energy shifts given by Eq. 11711 (dashed line) and Eq. USt 
(dashed-dotted line). The energy is measured in units of hui 
and the core width is L = 2<2ho- 



smooth potential, i.e. a spherical wave with wave vector 



feo = \/2[iE/h . This yields the relation 



8(E) = arctan 



— tan(fcL) 
k 



k () L. 



(16) 



In Fig. we show results for the energy of a few states 
with many nodes, obtained by direct calculation. The 
energy levels of our model are closely reproduced by Eq. 
©, while the predictions given in Ref. |||, 



ipout(r) oc e 



Vv 

Oho 



3 — r7 3 



4 

1 — 77 1 



2'< 



ho 



(15) 



where r\ — 2E/fko, \Fi(a,j 1 z) is the confluent hyperge- 
ometric function of the first kind 0] and the constant 
C = T (i-p) /2T (^-2) is chosen to ensure that ip out (r) 
vanishes at infinity |l6L Il7j . The allowed energies are 
found by matching at the boundary r = L the logarith- 
mic derivatives of the solutions inside and outside the 
well. 

The phase shift 6 for a given energy E may be found 
by equating at r = L the logarithmic derivatives of vpin 
and of the solution outside the core in the absence of the 



AE = - 



2-kTi 



tan 5 



\m\ 2 



(17) 



an 



d in Ref. [IJ, 



AE 



2irh z 
A* 



sin S cos 6 



(18) 



are correct only for small values of \S\ (modulo n). Note 
that if \S\ > ir/2 and \S/tt — n\ <^ 1 (where n is a suitable 
integer), the theoretical values of the energy given by 
expressions ifTTIl and iflSll agree with an energy eigenvalue 
of the system, but one with a different number of nodes 
inside the short-range potential. 
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III. A GENERALIZED GROSS PITAEVSKII 
EQUATION 

In the numerical example above, we considered a state 
with many nodes. For applications to Bose-Einstein con- 
densates, the relevant wave numbers are usually small, 
and therefore it is interesting to look at the opposite case 
of small but nonzero wave numbers. When the Wigner 
threshold condition k\a\ 1 is violated, the s-wave phase 
shift is no longer linear in the wave vector, and its en- 
ergy dependence can be written for potentials that fall 
off more rapidly than r~ 5 for large r as 



k cot S 



a 2 



(19) 



where r c is the effective range of the interaction |lj| ■ For 
small phase shifts 5 ~ tan<5 — tan 3 5/3 and we find: 

--=a(l-g 2 k 2 ) (20) 



where we have introduced 



ar c 



92 = T 

The energy shift is then given by 



AE 



2ixh 2 a 
(J* 



[i - 92 k 2 } iv(o)r 



(21) 



(22) 



It is interesting to note that in the case of hard spheres 
of diameter a the boundary condition on the surface of 
the sphere implies that 5 = —ka for all k and therefore 
g2 (and all higher terms in an expansion of 6 in powers of 
k) should vanish. Since from Eq. II19H r c = 2a/3 for the 
hard-sphere potential, one indeed finds gi = from Eq. 
(J2IJ- Our result is to be contrasted with the expression 
g-2 = a 2 — ar c /2 based on the approximation l(T8)l for the 
interaction energy fill. 

As shown in Ref. [11], the energy shift can be inserted 
into the energy functional to obtain a generalization of 
the GP equation for the condensate wave function that 
contains a derivative term in the interaction energy: 

*^*= "|^V 2 + l/(r) + C/o(|*| 2 +52V 2 |vI/| 2 ) 

(23) 

Contributions from higher partial waves may be included 
through the addition of higher derivative terms, as de- 
scribed for Fermi systems in Ref. 0. 

To test the validity of this prediction, we have per- 
formed numerical integrations of the relative motion for 
the two-body problem considered above, with a harmonic 
trapping potential and a short-range square well poten- 
tial, but for the lowest state rather than the excited states 
considered in Sec. II. The results are compared with those 
obtained from analytical approximations. Figure[2shows 
that the inclusion of the effective range correction with gi 
given by Eq. l|2"T|l greatly improves the prediction given 
by the simple scattering length approximation as soon as 
the condition k\a\ <C 1 is violated. 
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Figure 2: Energy shift for the ground state (n = 0) of the 
combined potential discussed in the text. The exact results 
(crosses) are compared with the energy shifts given by the 
formula AE = 2nh 2 f\tp(0)\ 2 / fi, where / = a (circles), / = 
— tan<5/fc (dotted line), / = — S/k (solid line), by the result 
that includes the effective range correction, / = a(l — g2k 2 ) 
(dashed-dotted line), and by the result of Ref. [Tl| . / = a(l — 
g-zk 2 ) (dashed line). The energy is measured in units of %lj 
and the core width is L = 0.25ah o . 



IV. APPLICATION TO RYDBERG 
MOLECULES 

In Rydberg atoms the valence electron is in a highly ex- 
cited state with principal quantum number n > 20, and 
within quantum defect theory it is described in terms of 
generalized hydrogenic wave functions. By analogy with 
the structures created around positive ions in liquid he- 
lium H} or BECs HHH, a Rydberg atom with a large 
electric dipole moment inside a BEC may create remark- 
able deformations of the condensate density in its sur- 
roundings. The additional degree of freedom introduced 
by the permanent electric dipole moment could also be 
used to realize conditional logic gates for quantum infor- 
mation processing j2^|. 

On the experimental side, Rydberg atoms have already 
been created in an ultracold environment [24], and there 
is a proposal to excite and trap single Rydberg atoms 
inside a BEC [25]. Being overall electrically neutral, 
Rydberg atoms are not accelerated by the stray elec- 
tric fields which are unavoidable in experimental vacuum 
chambers. This is particularly relevant since, in a typical 
apparatus for ultracold atoms, stray fields would drag an 
ion outside the condensate in a time much less than 1 
ms, making difficult the observation of the induced den- 
sity disturbances. 

In recent work it nas been predicted that the 

tailored excitation of single atoms in a BEC towards a 
Rydberg state would induce the formation of molecules 
characterized by ultra-long ranges (R ~ 2000 a.u.) and 
very large permanent electric dipole moments. In these 
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papers, the s-wave molecular potential between a ground 
state atom and a Rydberg atom was taken to be 

fi k 

where r and R are the positions of the electron and 
of the ground state atom relative to the Rydberg ion, 
5o[k(R)] is the energy-dependent phase shift and the elec- 
tron wave number k(R) is given by the hydrogenic rela- 
tion k 2 {R)/2 - l/R = -l/2n 2 . The authors of 0, E3 
follow Omont ja| and conclude that the appropriate po- 
tential for the excited electron is given, in the Born- 
Oppenheimer approximation, by 

U S (R) = E nl - g^ tanWfl)] ^ a 
/i k 

where E n \ and ip n io are the unperturbed atomic Rydberg 
energy and wave function (with quantum numbers n ~ 
30, I < 2 and m = 0). As discussed above, while the 
potential in Eq. l|2"3j) reproduces accurately eigenenergies 
and wave functions of the Schrodinger equation, it should 
not be used in mean-field calculations, where it yields 
unphysical divergent energy shifts in the presence of a 
resonance. 

This issue is particularly relevant for the scattering 
in the p-wave channel, where the e-Rb(5s) scattering 
phase shift 5i = i has a resonance at an energy of ap- 
proximately 30 meV, corresponding to a distance be- 
tween the Rb ion and the ground state atom of about 
700 a.u.: here Omont's expression for the energy shift 
AE(R) cx |V^„; (-R)| 2 tan<5i(fc)/fc 3 diverges. However, 
according to the arguments we have given, the appropri- 
ate effective interaction for a mean-field calculation is not 
given by this expression, but rather by one with the tan- 
gent of the phase shift replaced by the phase shift itself, 
and consequently there are no divergence problems. 

V. CONCLUSION 

In this paper we have studied the expression for the 
effective two-body interaction to be used in mean-field 



calculations of the energy of a state. All results agree 
in the limit of zero energy, but there are differences at 
higher energies. We have argued that the appropriate 
effective interaction is proportional to the phase shift, 
rather than other expressions that have been suggested, 
and we show that this holds analytically for the states 
of two particles in a trap in the quasiclassical approx- 
imation. We have calculated energy levels numerically 
for the problem of two particles interacting via a short- 
range square- well potential in the presence of a harmonic 
confining potential and have demonstrated that the ana- 
lytical expression in terms of the phase shift agrees well 
with the numerical data, both for the ground state and 
for excited states with many nodes. Since an effective 
interaction proportional to tan<5 gives the correct wave 
functions and energy eigenvalues when inserted in the 
Schrodinger equation, our results show that the choice of 
effective interaction depends on the application. 

In this article, we have also considered corrections to 
the Gross-Pitaevskii equation to allow for the nonzero 
energy of the relative motion of two particles and have de- 
rived a generalized Gross-Pitaevskii equation that takes 
into account the effective range. This equation gives a 
better approximation to the numerical results for energy 
eigenvalues than does an earlier proposal [ijj- A problem 
for future work is to include contributions from higher 
partial waves in the GP equation. Finally, we have ar- 
gued that the effective two-body interaction to be used 
in calculations of Rydberg atoms and molecules should 
be taken to be proportional to the phase shift. 
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